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Abstract 


Oja  (1983)  introduced  an  affine  invariant  multivariate  median.  Brown 
and  Hettmansperger  (1987)  developed  affine  invariant  rank  methods  in  the 
bivariate  location  model  based  on  Oja's  measure  of  scatter.  Analogs  of  the 
Wilcoxon  signed  rank  and  the  Mann-Whitney-Wilcoxon  rank  sum  tests  were 
presented.  In  .this  paper  we  discuss  scale,  rotation  and  affine  invariance. 

We~-bhen~drscuss  tests  in  the  bivariate  linear  model  for  an  overall  regression 

.  7„  ,  7  . 

effect/  In  particular  ws"  present  tests  in  the  one-  and  two-way  layouts  that 
are  analogs  of  the  Kruskal-Wallis  and  Friedman  tests. 


I 

! 

|  Key  words:  Affine  invariance,  AOV,  Kruskall-Wallis,  Friedman,  linear  models, 

spatial  ranks. 
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1.  Introd  -  Mon  and  Summary 


In  this  paper  we  present  affine  invariant  rank-like  teats  and  estimates 
in  the  bivariate  linear  model.  We  describe  a  general  test  for  an  overall 
regression  effect  and  then  specialize  the  test  to  the  one-  and  two-way 
layouts.  The  tests  are  affine  invariant  analogs  of  the  Kruskal-Wall is  (1952) 
test  and  the  Friedman  (1937)  test.  These  tests  are  based  on  a  notion  of  an 
affine  invariant  rank  vector  that  can  be  developed  from  Oja’s  (1983)  work  on 
an  affine  invariant  median.  The  analogs  of  the  Wileoxon  signed  rank  test  and 
the  Mann-Whitney-Wilcoxon  rank  sum  test  are  developed  by  the  authors  (1987). 

Along  the  way,  we  describe  three  possible  extensions  of  the  Lj  notion  to 
the  bivariate  setting.  The  three  types  of  bivariate  medians  that  result  are: 
the  scale  invariant  vector  of  component  medians,  the  rotation  invariant 
spatial  median,  and  the  affine  invariant  Oja  generalized  median.  These  form 
the  basis  for  the  extension  to  three  types  of  rank  vectors  in  the  linear 
model.  The  main  emphasis  will  be  on  the  affine  invariant  version  which  is 
illustrated  on  a  data  set  in  Section  3. 

We  indicate  in  Section  2  that  the  Oja  generalized  median  has  better 
efficiency  than  the  other  two  medians  for  an  underlying  normal  distribution. 
Although  we  do  not  have  asymptotic  efficiency  results  on  the  affine  invariant 
tests  in  the  linear  model,  we  illustrate  their  robustness  to  outliers  in  the 
example.  In  addition,  affine  invariance  may  be  an  important  consideration  in 
a  particular  application  and  the  tests  presented  here  provide  natural 
alternatives  to  the  traditional  tests  based  on  means. 
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2.  Bivariate  Medians  and  Invariance 

We  first  consider  the  development  of  rank  based  methods  in  the 
univariate  linear  model  from  the  Lj  norm  operating  on  a  single  sample.  For  a 
given  set  of  numbers  y ^ ,  ....  yn  define  their  center  by  0  such  that 

D(9)  =  Z  |y,-0|  ( 1 ) 

i  =  1  1 

is  a  minimum.  It  is  well  known  that  0  is  the  median  of  the  set  |y,,...,ynl, 
denoted  0  =  med  y  .  The  derivative  with  respect  to  0  exists  almost 
everywhere  and  is  given  by 

n 

Q ( 0 )  =  Z  sgn ( 0-y • ) .  (2) 

i  =  1  1 

The  estimate  0  =  med  y-  can  be  considered  the  solution  of  Q(0)  =  0  where 
is  interpreted  as  a  sign  change  for  n  odd. 

Next,  Let  r; {p)  -  y.  -  x-’p,  i  =  l,...n  represent  residuals  in  a  linear 

model,  where  x,'  =  (x-,,...,x-_)  and  P'  =  IP  ,...,#).  Then  an  estimate  of  P 

'*■  i  11  1  P  1  P 

is  determined  by  minimizing 

t  n 

IT(fi)  =  Z  D(r  Ap)  )  (3) 

i=l  1  ~ 

=  I  I  I  |  r ,  ( ## )  -  rdf?)  | 

i<J  J  ' 

n  N+  1 

=  Z  l Rank (  r ,  ( ##  >  )  -  '-jy-Mr.lP) 

i=l  1  '  L  1  * 
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=  i  Z  Q(r,  (*?)  )  r,  <*>  . 

The  second  equality  relates  the  estimate  of  if?  to  the  Lj  norm  and  the  third 
equality,  achieved  after  a  bit  of  algebra,  relates  the  estimate  of  fi  to  the 
ranks  of  the  residuals.  Jaeckel  (1972)  proposed  the  rank  version  of  D*lf?)  as 
a  criterion  for  estimating  0.  McKean  and  Hettmansperger  (1976)  proposed 
using  the  reduction  in  D  * ( /? )  when  passing  from  a  reduced  to  full  model  as  a 

criterion  for  testing  hypotheses  in  the  linear  model.  The  fourth  equality 
results  from  Q  (  r  ^  ( )  )  =  r.Jsgn(r1(^l  -  r  j  <  ^ )  )  =  2  ( Rank  (  r  1  ( f? )  )  -  (N+ll/2).  We 

consider  Q(r^(jf?)|  to  be  a  centered  rank  or  quantile  of  r^/?)  in  the  sample. 

Note  that  the  dispersion  measure  D*(|ff)  is  a  linear  function  of  the  residuals 

with  weights  determined  by  the  quantiles. 

Our  goal  is  to  carry  an  Lj  notion  over  to  the  bivariate  setting.  We 

then  will  have  a  bivariate  sample  median  with  efficiency  and  robustness 

properties  related  to  the  univariate  median.  Then  we  extend,  via  the  ideas 

in  (3),  to  the  linear  model.  The  gradient  Q(0)  acts  as  a  bivariate  quantile 

and  the  result  is  rank-like  tests  and  estimates  in  the  bivari ite  linear 

model.  The  bivariate  case  is  more  complex  than  the  univariate  case  and  we 

will  consider  three  different  Lj  notions. 

T 

let  =  <  y  i  1 , y  i  2 1 ’  1  "  I  »  ■  * • <  n  be  a  set  of  bivariate  observations,  and 
*  T 

define  v-  :  ( - y  i  £ » y  1 1 ) .  1  -  l,...,n.  We  now  present  three  bivariate 

dispersion  measures  each  of  which  entails  different  invariance  properties  for 
the  resulting  estimates  and  tests.  Let  9  -  I6| ,  then 


7VV  /  /  * V'V 
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n 


D  i  ( 0) 

=  i?l 1 iyil~®1 

.1  ♦  iyi2-02|} 

(  ’  ) 

n 

0  o 

1/2 

D2(0) 

=  i?1Uyil"°1>  +  (yi2'S2!  1 

(5) 

i'll  yJl  el  ‘ 

(6) 

D3  ( 0 ) 

=  2  2  | det 

i<j 

yi2  yj2  02 

1  • 

1  1  1 

The  respective  gradient  vectors  are  given  by 

Q^(0)  =  (2  sgnOj-y^) ,  2  sgn( 02-y i2 )  ) 

2  2 

Q2(0)  =  2  (lyire/  +  (yi2-02>  I  <®“Xi> 

'  yil  yji  el'  .  . 

Q3(0)  =  2  2  sgnldet  yi2  yj2  02  H^-j^) 

lJ  12  1 

The  dispersion  Dj(0!  in  (4)  with  Q ^ ( 0 >  in  (7)  defines  componentwise  Lj 
methods.  The  estimate  0  =  (med  y^j,  med  y^2),  *3  3ca^e  invariant  but  not 
rotation  invariant;  see  Bickel  (1964).  The  dispersion  D2(0)  in  (5)  with 
Q2(0)  in  (8)  defines  the  spatial  median  which  is  rotation  invariant  but  not 
scale  or  affine  invariant;  see  Brown  (1983).  Finally,  the  dispersion  D3(0) 
in  (6)  with  y3(0)  in  (9)  is  the  sum  of  areas  of  triangles  formed  by  taking 

and  pairs  of  observations  as  the  vertices.  Minimizing  the  sum  of  areas 
results  in  the  generalized  median  0  proposed  by  Oja  (1983).  The  Oja 

generalized  median  is  affine  invariant. 

Both  Dj(0)  and  D2<0)  reduce  to  the  Lj  norm  in  one  dimension.  This  is 


(7) 

(8) 

(9) 


f  CD 


not  the  case,  however,  with  ( © > •  It  seems  that  D ^ ( © )  is  more  intrinsically 

bivariate  than  the  others.  The  three  bivariate  medians:  componentwise 
medians,  spatial  median,  and  Oja  generalized  median,  have  somewhat  different 
efficiency  properties,  but  all  are  related  to  the  efficiency  of  the 
univariate  median.  The  univariate  median  has  efficiency  .637  relative  to  the 
mean  when  the  underlying  distribution  is  normal. 

Following  Bickel  (1964,  p . 1083) ,  we  define  the  efficiency  of  the 

bivariate  median  relative  to  the  bivariate  mean  to  be  [ g-var ( mean ) / 

1/2  222 

g-var (median ) ]  '  ,  where  g-variance  is  Og  (1-p  )  and  the  parameters 

are  taken  from  the  asymptotic  covariance  matrix.  Bickel’ s  Theorem  5.1 
provides  a  formula  for  the  efficiency  of  the  componentwise  medians  when 
the  underlying  distribution  is  bivariate  normal.  The  efficiency,  independent 
of  the  variances,  declines  as  a  function  of  the  correlation  coefficient  p  from 
.637  when  p  =  0.  Some  values  of  (p,  efficiency)  are  (0,  .637),  (.2,  .629), 
(.4,  .605),  (.6,  .558),  (.8,  .473)  and  (.9,  .396).  The  spatial  median,  on  the 
other  hand,  has  efficiency  that  does  depend  upon  the  variances  but  not  on  p. 
Using  Table  1  of  Brown  (1983)  and  Bickel’s  definition  of  efficiency,  it  is 
easy  to  see  that  the  spatial  median  is  generally  more  efficient  relative  to 
the  mean  for  an  underlying  bivariate  normal  distribution  than  the  vector  of 
componentwise  medians.  The  efficiency  depends  upon  the  ratio  of  standard 
deviations  of  the  two  components.  If  X  is  that  ratio,  then  some  values  of  (V, 
efficiency)  are:  (1,  .785),  (.8,  .783),  (.6,  .773),  (.4,  .747),  (.2,  .678), 
(.05,  .593)  and  (.01,  .321).  The  efficiency  deteriorates  as  the  contour 
ellipses  of  the  bivariate  normal  distribution  become  very  narrow.  Oja  and 
Niinimaa  (1985)  show  that  the  affine  invariant  Oja  generalized  median  has 
efficiency  relative  to  the  mean  equal  to  .785,  independent  of  the  variances 


and  correlation.  It  is  strictly  better  than  the  spatial  median  unless  the 


bivariate  normal  distribution  is  circular  in  which  case  the  efficiency  is  fhe 
same.  Hence,  the  Oja  generalized  median  has  superior  efficiency  properties  as 
well  as  enjoying  affine  invariance. 


3.  Tests  in  the  Linear  Model 


The  next  step  is  to  extend  these  Lj  notions  for  a  single  sample  to 
different  types  of  quantiles  in  the  linear  model.  Using  <7)  -  (9)  we  arrive 
at  the  vector  of  component  ranks,  the  spatial  ranks  and  the  at  fine  invariant 
ranks,  respectively.  These  ranks  then  provide  the  basis  for  tests. 


We  first  introduce  the  notation  of  the  bivariate  linear  model.  Let 


1  he 


an  n*2  observation  matrix  in  which  the  n  rows  are  independent  random  vectors 
such  that 


t  y  :  E 


T  1 


n 


Z  £, 


where  Z  is  an  n«(p+l)  matrix  of  known  regression  constants  and  is  a  '  p*  •. 


matrix  of  unknown  parameters. 

Let  r ,  ( ^  )  = 

Y 

-  l 

-  ?T£i 

denot e  the  i t  h 

residual  vector 

T 

where  z  is  the 

1 1  h  row  if 

V  ’ 

In  each  of 

the  three  cases 

(  1  )  -  i  t,  i  , 

i  t 

an  tie 

^hown  t  hat  ,  f ,  i 

t  " 

-  ZD  ir  i  I 
1  '  i  1  '  I  ' 
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1-52  D(A)i  >> 

dAst  '  A=0 


(15) 


where 


XTQ 


A  1  ' 
ql 


1  ) 


*i21 


,(  2 ) 


Tii  test  4  =  0  we  must  assess  the  size  of  L.  The  statistic  we  will 


jsi-  .s  'he  Hotel  1 ing-Lawley  trace  statistic  given  by 


S  :  (n-1  )Trace(LTIX^Xc)  1 L ( QTQ ) 1 


(16) 


-  i n-  1  i T race [ L ( QTQ )  1 LT ( xj Xc )  S 

T  T 

•  ne.o-  >.  is  the  mean  centered  design  matrix  and  we  assume  that  Q  Q  and  XcXc 

-ir-  r.  .risinguiar .  The  trace  form  of  S  is  not  immediately  intuitive  and  we 
»...  pr  o.  ide  two  motivations  for  S. 

: ti  'he  firs’  motivation,  we  roll  out  the  px2  matrix  L  by  columns  into 


,  i  1  I  ,  (  1  ),  (2)  ,  (2)  , 

1  •••Lp  L1  •'•Lp  >• 


(  17) 


: ;rs*  naif  of  ’he  vector  consists  in  the  statistics  relevant  to 
mpenent  and  similarly  for  the  second  half. 

•hat  ’ he  direct  or  Kronecker  product  of  two  matrices  is  defined 
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^r»s*Bt*u 


allB  *  ‘  •  alsB 


arlB 


ar  B 
rs 


rtxsu 


Now,  it  i9  straightforward  to  verify  that,  under  Hq :  A  =  0,  the  permutation 
covariance  matrix  of  LCQ^  is 


Cov(tcoi*  =  iPT  <9T9*'I&1 


Then  the  first  form  of  S  given  in  (16)  reduces  easily  to 


T"”  1  rn  1 

r\\  +  /  v  A  v 


s  =  ‘"-nLcoiU^)  *<?ac)  ]tcor 


(18) 


This  is  a  natural  quadratic  form  for  assessing  the  size  of  L.  Let  q1J  denote 

- 1  T  T 

the  (i,j)  element  of  (Q^Q)  and  let  )  from  (17),  then 

(18)  becomes 


S  =  (n-1)  I  i_qijL,i,T(xJxc) 


i=l  j=l 


T  -1 

Since  (Q  Q)  is  a  2*2  matrix  it  is  easy  to  find  the  inverse  and  get 


-xplicit  values  of  q1J.  If  we  let  q—  =  Z  q^*q^'  be  the  (i,j)  element  of 

r  1/2 

W  y )  and  let  :  d 1 2^ * q  1  lq22 *  '  '  t^ie  correlation  between  the  columns  of 


4,  t  hen 
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s  = 


-  (n-_Ll_  ,J_  ,(1) 


1-r 


12 


r  1  f 

?  qj 1  -  ’ -c-c 


T  ”  ^ 
<XiX,J  L 


(11  +  2-^ 


ql lq22  * 


L{  1 >T(Xjxc 


-1 


(2) 


♦  Li2)T(xJx  )  1LI2)). 
q22  '  *c*c 


I  19) 


This  shows  that  the  multivariate  test  statistic  is  composed  of  the 
corresponding  componentwise  test  statistics  along  with  a  cross  component 
statistic.  The  combining  weights  depend  on  the  variances  and  covariances  of 
the  quantiles. 

For  the  second  motivation,  we  roll  out  L  by  rows  into 


-row  l  l  l  i  p  p 


Under  Hn:  A  =  0,  the  permutation  covariance  matrix  of  L  is 


row 


C°V|W  =  A  (£5c”'9V 


Then  the  second  form  of  S  in  (16)  reduces  to 


s  =  'n-11  5rJ0wI,xJxcl'1*<QTQ)'1l!:roW 


■  L  (J 


another  natural  representation  of  S  as  a  quadratic  form  for  assessing  the  s.ze 
of  L . 

Let  denote  the  (i,j)  element  of  <  X c )  then  S  in  20  can  be 
written 

P  P  j  j  t  r  'I 

S  =  <n-l>  Z  Z  x'JL *  IW^l  L,  n 

1=1  j=l  c  1 1  '  '  '(Ji 

where  L  |  x  j  =  (L*-^L^*).  This  provides  a  nice  interpretation  because  •  n- 

square  of  the  univariate  statistic  is  replaced  by  a  quadratic  f  ra  ,n  ‘tie 
components  of  the  bivariate  statistic. 
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We  finally  note  that  to  teat  Hg:  A  =  0  vs  HA:  A  *  0,  that  is,  test  for 
an  overall  regression  effect,  we  reject  Hq  if 

S  >  t^(2p)  (22) 

O 

where  i^(2p)  is  a  chi-square  critical  value.  This  test  has  approximate 
significance  level  <*.  The  limiting  chi-square  distribution  follows  from 
standard  permutation  arguments;  see,  for  example,  Puri  and  Sen  (1985). 

We  next  specialize  the  statistic  S  to  the  one-  *'d  two-way  layouts.  In 

the  one-way  layout,  we  suppose  that  we  have  k  samples  of  size  n, . n^, 

from  continuous,  bivariate  distributions  with  location  vectors  0,,...,®^,  and 

we  wish  to  test  HQ:  0 1  :  ...  =  0^ .  The  design  matrix  X  consists  of  columns 

of  zeros  and  ones.  In  the  regression  setting  one  column  can  be  dropped  to 
insure  full  rank.  If  we  drop  the  first  column,  then  the  matrix  L,  rolled  out 

Dv  columns  is  :  i  L.!, '  ’  .  .  .  l!  1  ’  L ^  1  .  .  .  l!  2  1  )  and,  for  example,  L^1  is  the 

w(-  0 1  «,  K  L  K  L 

sum  of  the  first  components  of  the  combined  sample  quantiles  corresponding  to 

I  *  ( 

’he  second  sample.  The  matrix  Xr  along  with  (XcXri  is  given  by 

tie’  ’  mansperger  .S84,  p.258i.  The  element  x^ -  n  *  ♦  n,'  if  l  =  )  and  x^  1  = 
n  if  i  *  We  an  now  describe  *he  *wo  \ersions  of  S,  (  1 9 i  and  i 2 1  i  .  for 

•  he  ne-»ay  .ay  . ut  .  from  Is  we  f  .d 


H  H  H 

wnere  H  :  Z  1.  '  .  ft  -  Z  r.  I  and  rl  ‘  *  -  Z  n  ‘  ,  i  •  .  \'-tf 

'  ha'  n-  .  i  -  ’  .  s  i  nrustta.  -Wa  .  .  .  »  .^52  *  ■.  pe  «,» at  ;  sf  mpu'e.j  n  ♦  f,e 
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- 1  ( 2 )  (12) 

first  component  of  the  quantiles,  similarly  for  ( n— 1 ) q2 2  ^  *  "hile  H 

is  a  cross  component  statistic. 

From  (21)  we  find 

S  :  ( n— 1 )  Z  n"1  (l{1,lJ2))(QTQ)  1 (l{  1 >l| 2 } ) T.  (24) 

i=l  1  11  *  '  11 

This  corresponds  to  a  Kruskal-Wallis  calculation  using  the  k  quadratic  forms. 

If  q^*  =  -  (n+l)/2,  i  =  1,2,  j  =  l,...,n  are  the  centered 

componentwise  ranks  in  the  combined  sample,  then  S  is  the  test  discussed  by 
Puri  and  Sen  (1971,  p . 186 )  and  is  a  direct  generalization  of  the  univariate 
Kruskal-Wallis  test.  If  gj  is  taken  to  be  a  bivariate  spatial  quantile,  then 

a  rotation  invariant  test  results.  Our  primary  interest  is  in  describing  the 
affine  invariant  version  and  to  that  purpose  we  now  present  an  example. 
Example: 

In  Table  1  we  present  data  on  tail  length  and  wing  length  for  4 
subspecies  of  birds,  and  we  wish  to  carry  out  a  test  to  see  if  the  location 
vectors  differ  among  the  subspecies. 

-  Table  1  about  here  - 

Since  the  data  is  arrayed  in  a  one-way  layout,  we  will  use  the 
affine  invariant  test  based  on  S  in  (22).  We  find  S  =  36.5,  and  if  o  :  .01, 
the  approximate  critical  value  is  i'qj(6)  ^  16.8.  Hence,  the  null  hypothesis 
of  no  difference  among  the  4  population  location  vectors  is  easily  rejected. 

In  Figure  la,  we  present  plots  of  the  data  and  the  corresponding 
quantiles.  To  illustrate  the  effects  of  an  outlier,  we  changed  the  first 
component  of  the  first  observation  in  Table  1  from  207  to  2.7.  The  resulting 
value  of  S  is  33.0,  still  highly  significant.  In  Figure  lb,  we  present  the 
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new  plots  of  data  and  quantiles. 

-  Figure  1  about  here  - 

Generally,  the  statistic  S  will  be  less  effected  by  outlying  data  than 
the  normal  theory  based  Hotelling-Lawley  test.  The  authors  (1987)  provide  a 
two-sample  example  in  which  the  introduction  of  an  outlier  has  little  effect 
on  the  quantile  test  but  switches  the  traditional  test  from  significance  to 
nonsignificance . 

We  now  turn  to  the  two-way  layout.  We  consider  a  randomized  block 
design  and  describe  an  analog  to  the  Friedman  (1937)  test.  Suppose  we  have  k 
treatments  and  n  blocks.  In  the  univariate  two-way  layout  the  observations 
are  replaced  by  their  ranks  within  blocks.  Then  these  ranks  are  combined 
into  a  quadratic  form  based  on  the  k  treatment  rank  sums.  In  the  bivariate 
two-way  layout  we  replace  the  observation  vectors  by  the  bivariate  quantiles 
computed  within  blocks. 

Let  gTj  =  (qjj1,  dy1*  denote  the  quantile  in  the  (i,j)  cell, 
corresponding  to  block  l  and  treatment  j.  Then  =  (L!jl,Lji,l  = 

^  (  1  )  ^  I  2  ) 

(  Z  q,  :  Zqit)  is  the  sum  of  quant  lies  for  the  jth  treatment ,  j  =  1  ,  .  .  .  ,  k  . 
i=l  1J  i=l  1J 

It  remains  to  determine  the  permutation  covariance  matrix  of  the 
quantiles.  This  is  done  by  estimating  within  each  block  separately  then 
averaging  over  blocks.  We  define 

i  T  ,  n  k  r 

n (' k -  1  )  ^  ^  n ( k -  1  )  ,  - 1  I- i .j '  l“'M 


The  natural  analog  of  Friedman's  statistic  i. 
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S  =  (k-1 )  Z  ( l| 1 >  lJ  2 } ) <QTQ)  1(l[1)l|2))T  (26) 

i=l  1  1  -  -  ii 

k-1  .  1  F(l)  2q12  F ( 12 )  .  1  f(2), 

‘  7^7  ql  i  '  ‘  qHq22  F  q22  '  } 

where  F(j)  =  Z  ( L  J J 1 ) 2,  j  =  1,  2,  F(12)  =  Z  l[1,l}2),  and  qn  is  the  (i,j) 

1=1  1  i=l  1  1  1J 

T  1/2 

element  of  Q  Q  in  (24),  and  rjg  :  Q 1 2^ q  1  lq22 * 

If  ,  s  =  1,  2,  are  the  componentwise  centered  ranks  in  cell  (i,j), 
then  we  get  the  scale  invariant  test  described  by  Puri  and  Sen  (1971,  p.279). 
If  we  use  spatial  quantiles,  we  get  a  rotation  invariant  test  and  if  we  use 
the  affine  invariant  quantiles,  we  get  an  affine  invariant  test.  Standard 
permutation  theory  shows  that  an  approximate  size  a  test  rejects  the  null 
hypothesis  of  no  treatment  effect  if  S  >  x*(k-l). 


...it.**  < 


UC«J tldlc*  Im 


,v.  /.V^.v.n.^.v.v.v.V. 


■*.  t.  ■ 


References 


Bickel,  P.  J.  (1964).  On  some  alternative  estimates  of  shift  in  the 
P-variate  one  sample  problem,  Ann.  Math.  Statist..  35,  1079-1090. 

Brown,  B.  M.  (1983).  Statistical  uses  of  the  spatial  median.  J.  Rov. 
Statist.  Soc.  B,  45,  25-30. 

Brown,  B.  M.  and  Hettmansperger,  T.  P.  (1987).  Affine  invariant  rank  methods 
in  the  bivariate  location  model.  J.  Rov.  Statist.  Soc.  B,  to  appear. 

Friedman,  M.  (1937).  The  use  of  ranks  to  avoid  the  assumption  of  normality 
implicit  in  the  analysis  of  variance.  J.  Am.  Stat ■  Assoc.  32,  675-701. 

Hettmansperger,  T.  P.  (1984).  Statistical  Inference  Based  on  Ranks. 

Mew  York:  Wiley. 

Jaeekel,  L.  (1972).  Estimating  regression  coefficients  by  minimizing  the 
dispersion  of  the  residuals.  Ann.  Math.  Statist . .  43,  1449-1458. 

Kruskal,  W.  H.  and  Wallis,  W.  A.  (1952).  Use  of  ranks  in  one  criterion 
variance  analysis.  J.  Am.  Stat.  Assoc.  57,  583-621. 

McKean,  J.  W.  and  Hettmansperger,  T.  P.  (1976).  Tests  of  hypotheses  based  on 
ranks  in  the  general  linear  model.  ('omm.  Stat . -Theory  and  Methods 
A  5 i 8  )  ,  693-709. 

Oja,  H.  (i983).  Descriptive  statistics  for  multivariate  distributions. 

Stat.  and  Prob.  Letters,  l,  327-332. 


Oja,  H.  and  Siinimaa,  A.  (1985).  Asymptotic  properties  of  the  generalized 

median  in  the  case  of  multivariate  normal i t  >  .  J .  Roy .  Statist.  Soc ■  B, 

17,  372-377. 

Puri,  M.  1  .  and  Seri,  I' .  k.  (1971).  Nonuaranet  r  ic  Methods  in  Multivariate 
Vna lysis.  Sew  York :  Wiley. 

Pori,  M.  I  .  and  Sen,  1'.  k.  ■  isHYi.  Sunparamet  r  u Methods  in  denerai  linear 

Models.  John  Wilev,  \ev>  York. 


OBS  HIDDEN 


issBaaaafi 


Sj'Jjt" 


;  ->ivT  ACCci^iOM  ^O.j  i  «E:;piEht*; 


Invariant  Tests  in  Bivariate  Models  and 
the  Lj  Criterion 


<  Pfc.  '■jF  .»£PO«T  i  i  (.  Ml 


* .  3£fl-T"MihG  C  A  S  «£=."■ 


*  L  ■O  1 , 


B.  M.  Brovm,  University  of  Tasmania 
Thomas  P.  Hettmansperger ,  The  Pennsylvania 
State  University 


N00014-80-C-0741 


